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ABSTRACT 

After giving a summary of the basic-theoretical con- 
cept of quantization of the electromagnetic field in 
the presence of dispersing and absorbing (macro- 
scopic) bodies, their effect on spontaneous decay 
of an excited atom is studied. Various configura- 



tior s such as bulk material, planar half space me- (1987|); [John et al. (1994Q ; 



reabsorbed and reemitted. Obviously, the pho- 
tonic density of states can be modified by the pres- 
ence of macroscopic bodies, which, in the simplest 
case, change the boundary conditions for the elec- 
tromagnetic field. For the last years, engineering 
periodic dielectric structures (photonic crystals) has 
been of increasing interest flJohn (1987|); |Yablonovich| 

IKofman et al. (1994|); 



dia, spherical cavities, and microspheres are consid- 
ered. In particular, the influence of material absorp- 



Joannopoulous et al. (199q ); 



Soukoulis (Ed., 1996|) ; 



Woldcyohannes et al. (1999|); [Nikolopoulos et al\ 



tior j on the local-field correction, the decay rate, the 
line shift, and the emission pattern are examined. 
Further, the interplay between radiative losses and 
losses due to material absorption is analyzed. Fi- 
nally, the possibility of generating entangled states 
of two atoms coupled by a microsphere-assisted field 
is discussed. 



1 INTRODUCTION 

Spontaneous emission of an excited atom is not an 
immutable property of the atom, but it sensitively 
depends on the photonic spectral density of states 
that are involved in the atomic transition at the cho- 
sen location of the atom. Already Purcell (1946| ) 
pointed out that spontaneous emission can be en- 
hanced when the atom is inside a cavity and its 



(2000 ); [Zhu et al. (20001 ); [Schriemer et al. (20011 )].' 



It is worth noting that spontaneous emission may 
be regarded as being a basic process in the rapidly 
growing field of cavity quantum electrodynamics 
(QED), where strong (resonant) interactions of a sin- 
gle or a few atoms with a single or a few radiation- 
field modes formed by material bodies are stud- 
ied. Cavity QED itself has offered novel possibilities 
of testing fundamental aspects of quantum physics, 
such as quantum nondemolition measurement, com- 
plementarity, and entanglement [for reviews, see, 
e.g., [Hinds (1991]); [Haroche (1992Q; [ Meschede (19921 ); 
Meystre (1992|) ; |Berman(Ed., 1994 ); paroche(1998[) ; 
Kimble (199$ ; [Walther (199^) 1- 

Spontaneous emission in the presence of material 
bodies is not only interesting from the point of view 
of fundamental research, but it has also offered a 
number of interesting applications. It can provide a 
reliable and efficient single-photon source to be used 



transition is in resonance with a cavity mode. The 



opposite case of inhibition of spontaneous emission 
is also possible | Kleppner (1981| )1- It is further well 



in quantum information processing | De Martini et 
al. (1996j ), |Kitson et al. (1998| )]. The sensitivity to 



the ambient medium of resonance fluorescence is cru- 
cial in scanning near-field optical microscopy [ Betzig 



kn own that the decay process can even be reversed et aL (1991) ; |Kopelman et al. (1993D ; |Bian et al. 



by strongly coupling the atom to a 
cavity-field mode, so that the emitted photon can be 



sumcientiy snarp (1995| ); [Henkel et al. (199% |Gersen et al. (20001 )]. 



Another important potential application is the so- 



called thrcsholdless laser [Dc Martini et al. 



'On leave from the Institute of Physics, National Center 
for Sciences and Technology, 1 Mac Dinh Chi Street, District 
1, Ho Chi Minh City, Vietnam. 



Yamamoto et al. (1993); Protsenko et al. (1999)]. 
In a conventional laser, only a small portion of the 
spontaneous emission is channeled into the lasing 
mode formed by the cavity mirrors, the rest being 



1 



lost to the free space modes. In a microcavity, due 
to strongly modified emission pattern and enhanced 
emission rate, a large portion of spontaneously emit- 
ted photons is stored in the cavity resonance mode. 
Losses due to excitation of free-space modes are thus 
drastically reduced and ultralow threshold lasing can 
be achieved. 

Controlling of the spontaneous decay also plays 
an important role in solid-state systems, where in- 
stead of atoms, quantum well or quantum dot exci- 



tons play the role of the emitters [Ybkoyama et al 



(Ed |s., 1995| ); [Khitrova et al. (1999| ); [Yamamoto"d 
al. (20001) ]. So, the improved directionality of the 



spontaneously emitted light could have dramatic im- 
pact on manufacture of high-efficient light-emitting 
diodes and displays | Yamamoto et al. (1993| )], and 
the spectral narrowing could help to increase the 
transmission capacity of optical fiber systems where 



chromatic dispersion is the limiting factor [Hunt et 
al 



Although certain properties of spontaneous emis- 
sion such as the decay rate can be described clas- 
sically, using the model of a classically oscillat- 
ing dipole interacting with its own radiation field 
[|Chance et al. (1978Q ; |Wylie et al. (1985Q ; [E 
och 3(1992| )], spontaneous emission is an intrinsically 
quantum mechanical process. Its proper description 
requires quantization of both the atom and the ra- 
diation field. Obviously, in the presence of mate- 
rial bodies the medium-assisted electromagnetic field 
must be quantized. As long as the medium can be 
regarded as being nondispersing and nonabsorbing, 
whose (real) permittivity changes with space in gen- 



eral, eiectromagnetic-Ueld quantization can be per- 



formed, using, e.g., generalized orthogonal-mode ex- 
| Knoll et al. (1992]) ]. However, the concept 



pansion 



material arises. Obviously, spontaneous decay must 
not necessarily be accompanied by a really observ- 
able photon, if the atom is near an absorbing body, 
and the question is of what is the (average) fraction 
of emitted light. Another question is of how can 
absorption modify the local field felt by the atom. 
A rigorous approach to the problem has acquired 
even more significance with the recent progress in 
designing certain types of microcavities (e.g., mi- 
crospheres), where the ultimate quality level deter- 
mined by intrinsic material losses has been achieved 
[ Porodetsky et al. (1996| )]. 



Although there has been a large body of theo- 
retical work on medium-assisted spontaneous emis- 
sion, material absorption has been ignored usually. 
Roughly speaking, there have been two concepts to 
treat absorption, namely the microscopic and the 
macroscopic approach. The microscopic approach 
starts from a microscopic model o f the medium [[Leej 
et al. (1995| ); [Yeung et al. (199CQ ; |Juzeliunas (1997D ; 
Flcischhaucr (19991 ); |Crenshaw et al. (2000"i| |b|); 
Wubs et al. (2001 )]. Accordingly, the underlying 



total Hamiltonian typically consists of the Hamil- 
tonians of the free atom, the free radiation field, 
the atomic systems of the medium, and all the mu- 
tual interactions. The resulting equations of motion 
for the coupled system are then tried to rewrite in 
order to eliminate, on applying various approxima- 
tion schemes, the medium variables and to obtain 
closed equations of motion for the atom-field sys- 
tem only. In this way, the life time of an excited 
atom in absorbing bulk material Lee et al. (1995); 
luzeliunas (1997| ); |Fleischhauer (1999] ); prenshaw e\ 



al. (2000a| Jb|); |Wubs et al. (2001| )1, the initial tran- 
sient regime [ |Wubs et al. (2001|) ], and the problem 



fail* when material absorption is included and the 



(spatially varying) permittivity becomes a complex a l 



of local field corrections | Juzeliunas (1997 ); Fleiscb 
hauer (1999p ; |Crenshaw et al. (20004 p|); |Wubs et 
(2001 )] have been studied, and the problem 



of spontaneous emission of an excited atom near 



function of frequency. The systematic study of the 

pro blem during the last years has generalized earlier an absorbing interface has been considered [|Yeung 
res ults and onered poweriul methods to deal with et a j (1996p ]. The concepts typically borrow, at 



the special requirements of quantizing the electro- 
magnetic field in absorbing media [for a review, see 



Knoll et al. (2001 )]. 

There are many reasons why inclusion of material 
absorption in the study of spontaneous emission is 
desired. One might ask what would happen when 
the atomic transition frequency becomes close to a 
medium resonance, where absorption is strong. In 
particular the question of the effect of absorption on 
spontaneous emission in the presence of band-gap 



some stage of calculation, from macroscopic electro- 
dynamics, e.g., when a (model-specific) permittivity 
is introduced, boundary conditions at surfaces of dis- 
continuity are set or local-field corrections within the 
framework of cavity models are considered. Apart 
from the fact that the (simplified) microscopic mod- 
els do not yield, in general, the exact permittivities, 
the calculations can become rather involved, partic- 
ularly when surfaces of discontinuity are taken into 
account | Yeung et al. (1996j )]. Further, the elimi- 



2 



nation of the medium variables must be done very 
carefully in order to ensure that the equal-time com- 
mutation relations are preserved. If this is not the 
case | Crenshaw et al. (2000a , b|)l, the results are ques- 
tionable. 

In the macroscopic approach, the medium is de- 
scribed, from the very beginning, in terms of a spa- 
tially varying permittivity, which is a complex func- 



Hence, knowing the Green tensor of the classical 
problem for given complex permittivity, the decay 
rate and the line shift are known as well. Equations 



( |1.3[) and (|1.4j) were used in order to calculate decay 
rates and line shifts for an excited atom near a real- 



istic (i.e., absorbing) metallic sphere |Ruppin (1982|) ; 
Agarwal et al. (1983|) ], near an absorbing interface 
[|Agarwal (1975^; [Agarwal et al. (1977| ); |Wylie et at 



tior . of frequency, e(r, u>), that satisfies the Kramers- 
Kronig relations. This approach has - similar to 
classical optics - the benefit of being universally 
valid, because it uses only general physical proper- 



(1984 , 1985)], and in a planar cavity filled with an 



absorbing medium [Tomas et al. (1999)]. Based on 
Eq. (|1.3|) , spontaneous emission of an excited atom 
in absorbing bulk material was studied [Barnett et 



ties , without the need of involved ab initio calcu- al. (1992 )], including local field corrections jBarnett 



lati bns. Clearly, this concept is valid only to some et al. (1996| )1 . The associated line shift (without 

local-field correction) was considered in Welton's in- 
terpretation [[Matloob (200q )l. Spont aneous decay 
of an atom at the center of an absorbing sphere has 
been calculated in Tomas (2001|) (with local field 
correction) . 

In what follows we restrict our attention to the 
macroscopic approach that is based on QED in dis- 
persing and absorbing media, within the framework 
of a source-quantity representation of the electro- 
magnetic field in terms of the (classical) Green ten- 
sor of the macroscopic Maxwell equations and ap 
propriately chosen fundamental bosonic fields [Hut 



approximately fixed length scale which exceeds the 
average distance of two atoms. With regard to the 
calculation of the lifetimes and line shifts, the macro- 
scopic approach is simple. It is well known that, 
according to Fermi's golden rule, the rate of spon- 
taneous decay V of an excited atom [position r^, 
(real) transition dipole moment d, transition fre- 
quency loa\ can be expressed in terms of an electric- 
field correlation function as follows [see, e.g., 
(19©]: 



Loudon 



2tt 



[cf. 



[Ab 



d^(0 | dE(rA.^)(SET(rA.^)d | n) 



tner et al. (1992| ); |Ho et al. (1993[) ; |Matloob etal. 



Eu,5. (|2.l| ) (|Q.3|) ]. It is alau well knuwii [mt, g.g., 



ikosov el al. (1975] )] that, in agreement with Lhe 



TT991); |Gruner et al. (1996|); [Ho et al. (199Sj ); ^chiel 
et al. (1998| ); [Knoll et al. (200 1| )]. The quantiza- 



dissipation-fluctuation theorem, the relation 

(0|i;(r ) w)®^ t (r / J w')|0) 
fiuj 2 



7reoc z 



ImG(r,r',w)(5(w-w') (1.2) 



(c, vacuum velocity of light) is valid, where 
G(r, r' , to) is the Green tensor of the classical, macro- 
scopic Maxwell equations. Combining Eqs. ( |l.l| ) and 



(1.2) yields 



r _ M 



dImG(r A ,r A , wa) d 



(1-3) 



[see also [Agarwal (19751) ; jWylie et al. (198j [19850 ]. 
The line shift can be calculated in a similar way to 
obtain 



tion scheme is outlined in Section ||. In Section ||, 
the basic formulas for studying the spontaneous de- 
cay of an excited atom are given, which cover both 
the strong- and the weak-coupling regime. Section 
|] is devoted to the spontaneous decay in bulk mate- 
rial, with special emphasis on local-field effects. The 
problem of spontaneous decay near a planar inter- 
face is considered in Section S and Sections || and [5], 
respectively, present results for an atom in a spher- 
ical cavity and near a microsphere. In Section || 
a system of two atoms coupled to a microsphere is 
analyzed, with special emphasis on entangled-state 
preparation. Finally, a summary is given in Section 



2 QUANTIZATION SCHEME 



Sua 



V 



nhe, 



du 



u 2 dim G{ya, i\4, u) d 



oJo 



U — OJA 



[1A) 



[V, principal value; for a more rigorous derivation of 
Eqs. dL3|) and Q, see Section |3T2TT . 



Following [Ho et al. (199^ ); [Scheel et al. (1998|) ; 
Knoll et al. (2001), we first consider the electromag- 
netic field in the presence of dispersing and absorb- 
ing macroscopic bodies in the case where no addi- 
tional atomic sources are present. The electric-field 
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operator E can be represented in the form of 

E(r) = EW(r)+EH(r), (2.1) 
E(~)(r) = [E( + )(r)] f , (2.2) 

/>oo 

E (+) (r) = / dwE(r,w), (2.3) 
Jo 

and the induction-field operator B accordingly. The 
fields E and B then satisfy the macroscopic Maxwell 
equations 

VB(r,w) = 0, (2.4) 
e e(r, w)E(r, u) = p n (r,uj), (2.5) 



VxB(r,w)= W j N (r,w 



V x E(r,Lj) = iu)B(r, uj), 



(2.6) 

2 £ (r,u/)E(r,a/). (2.7) 



As already mentioned, the real part er and the imag- 
inary part Ei of the complex (relative) permittivity 
e(r, uj) satisfy (for any r) the Kramers-Kronig rela- 
tions. The operator noise charge and current den- 
sities p N (r,u;) and j (r, w) respectively, which are 
associated with absorption, are related to the oper- 
ator noise polarization P N (r,u) as 



J N (r,u;) 



-VP N (r,u;). 
-iujP N (r,uj), 



where 



P N (r.^) i\l^ei(r,u)f(r,u 



(2.8) 
(2.9) 

(2.10) 



Here, f(r, lS) and V(r,uj) are bosonic fields which 
play the role of the fundamental variables of the 
composed system (electromagnetic field and medium 
including a dissipative system), 

[f i (T,u),fj(T f ,u , )]=6 ij S(T-i f )S(u-u/), (2.11) 

[f i (r,Lj)J j (r',u')]=0. (2.12) 

From Eqs. (U) " (EDI) ^ follows that E can be 
written in the form 



E(r,w) = i\ 5" 

V vre c z 

x ! d 3 rVe/ (r', w) G(r, r', u)f (r', u), (2.13) 

and B = (iuj)^ 1 V x E accordingly, where G(r, r', u) 
is the classical Green tensor satisfying the equation 

,2 



s(t,u)) — V x Vx 



G(r,r',u) = -6{r-r') (2.14) 



together with the boundary condition at infinity 
[S(r), dyadic <5-function] . In this way, the electric 
field and the induction field are expressed in terms 
of a continuum set of the bosonic fields f(r, u) [and 
f'(r, u)], and the Hamiltonian of the composed sys- 
tem reads (without the infinite ground-state energy) 



H= /d 3 r / &uftu>$(T,u)f(T, 



(2.15) 



Using Eq. fl2.13| ) [together with Eqs. ( |2.1| ) and 
Q2.3[) 1, one can introduce scalar and vector potentials 
<p and A, respectively, and express them in terms of 
the fundamental bosonic fields. In particular, in the 
Coulomb gauge one obtains 

-V0(r)=E«(r), (2.16) 

f'°° - i 

A(r) = / dw(iw)- 1 E (r,w)+H.c, (2.17) 
Jo 

where 

E^ir) = J dV^l^r-rOEtr'), (2.18) 

with 5 _L (r) and <5"(r) being the dyadic transverse 
and longitudinal 5-functions, respectively. 

We now consider the interaction of the medium- 
assisted electromagnetic field with additional point 
charges q a . Applying the minimal-coupling scheme, 
we may write the complete Hamiltonian in the form 
of 

H = / d 3 r / duo Huj f t( r , w)f (r, w) 



q a A(r a ] 



E 2 



+ \ J d 3 rp A (r)^(r) + J &h p A {r)${r), (2.19) 

where r a is the position operator and p a is the 
canonical momentum operator of the ath charged 
particle of mass m a . The Hamiltonian ( 2.19| ) con- 
sists of four terms. The first term is the energy ob- 



served when the particles are absent [cf. Eq. ( 2.15 )]. 
The second term is the kinetic energy of the par- 
ticles, and the third and fourth terms are their 
Coulomb energies, where the potential (pA can be 
given by 



<PA(r) 



dV 



PA<J?) 



47T£o|r — r'| 



with 



(2.20) 



(2-21) 
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being the charge density. Obviously, the last term in 



Eq. ( 2.19| ) is the Coulomb energy of interaction of the 
particles with the medium. Note that all terms are 
expressed in terms of the dynamical variables f (r, u>), 
fT(r, lo), r a , and p Q . It is worth noting the quan- 
tization scheme is fully equivalent to the so-called 
auxiliary-field scheme introduced by Tip (1998| ) 
[for details, see Tip et al. (2001| )]. 



3 SPONTANEOUS DECAY: 
GENERAL FORMALISM 

3.1 Basic equations 

Let us consider N two- level atoms [positions ta, 
transition frequencies <jJa (A = 1, 2, N)] that res- 
onantly interact with radiation via electric-dipole 
transition (dipole moments (1a)- Let us further as- 
sume that the atoms are sufficiently far from each 
other, so that the interatom Coulomb interaction 
can be ignored. In this case, the electric-dipole ap- 
proximation and the rotating wave approximation 
apply, and the minimal-coupling Hamiltonian takes 
the form of flHo et al. (200C|) ; |Knoll et al. T200T|) 1 



H 



d 3 r 



du) fkv f^(r, uj)¥ (r, to) 



o 



^[^E(+)(r A )d A + H.c.]. (3.1) 

A A 
Here and in the following, the two-level atoms are 
described in terms of the Pauli operators a a, o\, 
and a Az . 

For a single-quantum excitation of the system, the 
system wave function at time t can be written as 

\m) = T, c u A ^~ i{UA ~ Q)t \ u A)\m 



+ 



A 

d 3 r 



x e 



duj 

o 



C L i{r,u},t) 
\L)\{Ur,co)}) 



(3.2) 



(fi = 2^A u a)- Here, \Ua) is the excited atomic 
state, where the ^4th atom is in the upper state and 
all the other atoms are in the lower state, and \L) 
is the atomic state, where all atoms are in the lower 
state. Accordingly, |{0}) is the vacuum state of the 
rest of the system, and |{lj(r,u;)}) is the state, where 
it is excited in a single-quantum Fock state. The 
Schrodinger equation yields 



yJei{Y,u) d A G(r A ,r,uj)C L (r,uj,t) j, (3.3) 



x ^ei(r,u)d A G*(r A ,r,uj)C UA (t). (3.4) 
We now substitute the result of formal integration 



of Eq. (|3.4[) into Eq. (|3.3|) . Making use of the rela- 
tionship 



Im%(r,r',w) = J d 3 



-n£l{S,V) 



x G ik (r,s,u)G* k (r',s,uj) 



(3.5) 



we obtain the following (closed) system of integro- 
differential equations: 

Cu A (t) = V f dt' K AA > (t, t') C Ua , (0 
Jo 



1 



^/ire hj 



du / d 3 r 



-i(itJ— U>A)t 



V 7 £i(r, u) d A G(r A , r, u)C L (r, to, 0) }, (3.6) 



K AA ,(t,t') 



■1 



du 



-i(u) — u)A)t 



x e 



TllTEQ j 



Note that 



K AA '(t,t') = K%, A (t',t), 



(3.7) 



(3.* 



because of the reciprocity theorem. 

The excitation can initially reside in either an 
atom or the medium-assisted electromagnetic field. 
The latter case, i.e., Cr,(r, u, 0) ^0 in Eq. (|3.6[) , 
could be realized, for example, by coupling the field 
first to an excited atom D in a time interval At such 
that, according to Eq. (|3.4[), Ci,(r, CJ, 0) reads 



C*(r >W) 0)=/° L= ^ ^ D)t ' 

J -At V7T£oA C z 

x y/ei(T,u)d D G*(T D ,r,u)Cv D tf), (3.9) 

where Cu D (t) describes the single-atom decay ac- 
cording to Eq. ( |3.13 ) given below. Substitution of 
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the expression (3.£) into Eq. (3.C) then yields 



CU A (t) = J2 I dt ' K AA> (t, t') C Ua , (t>) 
A' Jo 

+ / dt'K AD (t,t')C UD (t'). (3.10) 
J-At 



3.2 Single-Atom Decay 
3.2.1 Atomic Dynamics 

Let us restrict our attention to the spontaneous 
decay of a single atom. For the initial condi- 
tion Cl(t, cj, 0) =0, Eq. ( |3.6| ) becomes {Cy = Cu A , 
d = d A ) 

Cu(t)= f ' dt'K(t-t')Cu(t'), (3.11) 
J o 

where 

K(t-t')=K AA (t,t'). (3.12) 

Integrating both sides of Eq. ( |3.11[ ) with respect to 
time, we easily derive, on changing the order of inte- 
grations on the right-hand side, a Volterra integral 
equation of the second kind, 

Cu(t)= f dt' K(t-t')Cu(t') + l (3.13) 



[Cu(0) = 1], where, according to Eqs. fl3.12|) and 



K{t-t') 



■1 



doj 



1- 



-i(w—u>A)(t—t') 



TmeoJo 
uj 2 dim G(i\4, r^, w)d 



(3.14) 



It is worth noting that Eqs. ( 3.1 1| ) and ( 3.13| ) ap- 
ply to the spontaneous decay of an atom in the 
presence of an arbitrary configuration of dispersing 
and absorbing macroscopic bodies. All the matter 
parameters that are relevant for the atomic evo- 
lution are contained, via the Green tensor, in the 
kernel functions ( |3,12j ) and (3.14). In particular 
when absorption is disregarded and the permittivity 
is regarded as being a real, frequency-independent 
quantity (which of course can change with space), 
then the formalism yields the results of standard 
mode decomposition, obtained by Laplace transform 
techniques [ Lewenstein et al. (1988a , bh ; John et 
al. (1994)] and delay-differential-equation techniques 
flCook et al. (1987|) ; |Feng et al. (1989|) ; |Ho et al 



(1999j )]. It should be pointed out that the Green ten- 



sor has been available for a large variety of configura- 
tions such as planarly, spherically, and cylindrically 
multilayered media jTai (1994D ; |Chew (1995|) ]. 

In order to study the case where the atom is sur- 
rounded by matter, the atom should be assumed to 
be localized in some small free-space region, so that 
the Green tensor at the position of the atom reads 



G(r A , r A , uj) = G v (r A , r A , u) + G H (r A , r A ,tn 



Ri 



where G v is the vacuum Green tensor with 
-.V, - u 



lmG (r A ,r a ,uj) 



6lTC 



(3.15) 



(3.16) 



(Appendix |A.l| ), and G R describes the effects of re- 
flections at the (surface of discontinuity of the) sur- 
rounding medium. The contribution of G v to K in 
Eq. (|3.14[) can be treated in the Markov approxima- 
tion (see below), thus 



K{t - f) 



^0 + 



da; 



o Jo 



dim G (r A , r A , u>)d 
i(u - LO A ) 



-i(u>—u)j\)(t—t') 



(3.17) 



where Tq is the well-known decay rate in free space, 



rv 



u A d 2 

3hTT£oC 3 



(3.18) 



The integro-differential equation ( 3.11| ) [or the inte- 
gral equation ( |3.13j )] together with the kernel func- 
tion ( |3.12| ) [or (|3.17| )1 can be regarded as the basic 
equation for studying the influence of an arbitrary 
configuration of dispersing and absorbing matter on 
the spontaneous decay of an excited atom. 

Weak Coupling When the Markov approxima- 
tion applies, i.e., when in a coarse-grained descrip- 
tion of the atomic motion memory effects are disre- 
garded, then we may let 



((uj a - u) (3.19) 



i(uj A - oj) 

in Eq. (|3.14| ) [((x) = ir5(x) + iV /x], and thus 

K(t-t') = -\Y + i8u) A , (3.20) 



where T and 5uj a are respectively given by Eqs. ( |P 
and ( |1.4|) . Substitution of the expression ( [OOT) into 
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Eq. ( |3.13| ) for the kernel function yields the familiar 
(weak-coupling) result that 



Cu(t) = exp[(-±r + i5uj A ) t] . 



(3.21) 



Application of Eqs. (HI) , (HI) , (|3l5|) , and (|1|) 
yields 

2a; 2 

r = r o + ^2 dImG!/? ( r ^' r ^'^) d ' ( 3 - 22 ) 



heoc 2 



7T 



dReG R (r A ,r A ,LO A )d 
R 



00 n iv 2 dlmG R (r A ,r A ,uj)d 
da; — = — 



UJ + U A 



(3.23) 



In Eq. (3.23), the Kramers-Kronig relations have 
been used and the divergent contribution of the vac- 
uum to the line shift is thought of as being included 
in the atomic transition frequency oo A . It is not diffi- 
cult to see that in Eq. ( |3.23 ) the second term, which 
is only weakly sensitive to the atomic transition fre- 
quency, is small compared to the first one and can 
therefore be neglected in general. 

Strong Coupling When the atomic transition 
frequency approaches a resonance frequency of a 
resonator-like equipment of macroscopic bodies, 
then the strength of the coupling between the atom 
and the electromagnetic field can increase to such an 
extent that the weak-coupling approximation fails 



and the integral equation ( 3.11 ) must be consid- 
ered. Let us assume, for simplicity, that only a single 
(field-)resonance line of Lorentzian shape is involved 
in the atom-field interaction. In this case, the kernel 
function ( |3.12| ) may be approximated by 



K(t - t') 



r c (Aa; c ) 2 i{uJc -u A) (t-t>) 



2tt 



-i(u>-L> C )(t-t') 



duj ■ 



00 (uj — uJc) 2 + (AlucY 



-i(u c -u>A)(t-t')„-Au,c\t-t'\ 



(3.24) 



and thus the integral equation ( 3.1l|) corresponds to 
the differential equation [ [Ho et al. (2000| )] 

C v {t) + [i(u c - uj a ) + Au c ] C v (t) 

+ \YcAucCu{t) = 0. (3.25) 

Here, u>c and Aujc are respectively the mid- 
frequency and the line width of the field resonance 



associated with the bodies, and Tc is the (weak- 
coupling) decay rate at ujq- 

Equation ( [3.25 ) typically applies to the strong- 
coupling regime for an arbitrary resonator configura- 
tion, provided that the field that effectively interacts 
with the atom can be regarded as being a single- 
resonance field of Lorentzian shape.Q In particu- 
lar, when material absorption is disregarded, then 
the line broadening solely results from the radiative 
losses due to the input-o utput coupling JCook et al. 



fl987D ; [Lai et al. (198Sp ; [Feng et al. [jggg j] 



Equation ( 3.25| ) reveals that the upper-state prob- 
ability amplitude of the atom obeys the equation 
of motion for a damped harmonic oscillator. In 
the strong-coupling regime, where uj a = u>c and 
Q ^ Aljc, damped Rabi oscillations are observed: 



Cu(t) 



-Aui c t/2 



cos(fit/2), 



where the Rabi frequency reads 



O = y / 2r c A 



uj c ■ 



(3.26) 



(3.27) 



3.2.2 Emitted-Light Intensity 

It is well known that the intensity of light registered 
by a point-like photodetector at position r and time 
t is given by 

I(r,t) = <V^)|E(-)(r)E^)(r)|V^)>. (3.28) 

The emitted-light intensity associated with the spon- 
taneous decay of an excited atom in the presence of 
dispersing and absorbing matter can be obtained by 
combining Eqs. Q - @. @), ©, and (p|) . 
The result is 



J(r,t) 



7T£oC 



dt' / duj 



Cu(t') 



e -i( W -uu)(t-f ^ w j d 



(3.29) 



where, in the spirit of the rotating wave approxima- 
tion used, J 1 = uj\ has been set in the frequency inte- 
gral. Again, all relevant matter parameters are con- 
tained in the Green tensor. In contrast to Eq. ( 3.13| ) 
[together with the kernel function (|3.17]) 1, Eq. ( 3.29Q 
requires information about the Green tensor at dif- 
ferent space points. In particular, its dependence 



1 Equations of the type ( p.25| ) can also be obtained within 
the framework of standard (Markovian) quantum noise theory, 
where an appropriately chosen unda mped mode is c oupled to 
a two-level atom and some reservoir |Sachdev (1984)]. 



7 



on space and frequency essentially determines the 
retardation effects. 

In the simplest case of free space we have 



lmG V (r,r A ,uj) 

x I e luJp/c - 



1 



8inp \ p" 
-™p/A + 0(p~ 2 ) (3.30) 



(p = y — y a ; Appendix |A.1|) . We substitute 
Eqs. ( ^2l|) (r = r ) and (Q) into Eq. (|3T29| ), cal- 
culate the time integral, and extend the lower limit 
in the frequency integral to — oo, 



du \J"ph- e -i"p/c\ 



-(T Q /2+iuj' A )t_ -iut 



i[u-(u' A -iT /2)} 

P 



= -2vrexp {-\Y Q -iu' A ) (t- ^ e(t 
[@(x), unit step function], where 

uJ A = °JA — 8lj a ■ (3.32) 
Thus, the well-known (far- field) result 



I(r,t) 



u A ds'm6 
4tteoc 2 p 



-r (t- P /c) ®(t-p/ c ) (3.33) 



is recognized (8, angle between p and d). 

It should be noted that the general expression 
( 3.29| ) is valid for an arbitrary coupling regime. In 
particular, in the weak-coupling regime the Markov 
approximation applies, and Cu{t') can be taken 
at t' = t and put in front of the time integral in 
Eq. ( 3.29| ), with Cu(t) being simply the exponential 
( |3 .21 ). Equation ( |3.29| ) thus simplifies to 



I(r,t) ~ |F(r,r A ,u; A )| 2 e 



-rt 



(3.34) 



where 



F(y,y a ,uj a ) ~- 



IUJ, 



G(Y,Y A ,U} A )d 



du 



£qC 

Im G(r, r A , uj)d 



(3.35) 



7T Jq U + LO A 

Since the second term on the right-hand side of 



Eq. (3.35) is small compared to the first one, it 
can be omitted, and the spatial distribution of the 
emitted-light intensity (emission pattern) can be 
given by, on disregarding transit time delay, 



\F(r,r A ,u A ) 



e c z 



G(Y,Y A ,UJ A )d 



(3.36) 



Material absorption gives rise to nonradiative de- 
cay. The fraction of really emitted radiation energy 
can be obtained by integration of the Pointing- vector 
expectation value with respect to time and over the 
surface of a sphere whose radius r is much larger 
than the extension of the system consisting of the 
macroscopic bodies and the atom, 

/'OO f2ir rn 

W = 2ce dt d(j) d6r 2 sm9 1(r,t). (3.37) 



o 











The ratio W/Wq {Wq = hu> A ) then gives us a mea- 
sure of the emitted radiation energy, and accord- 
ingly, 1 — W/Wq measures the energy absorbed by 
the bodies. 



(3.31) 3.2.3 Emitted-Light Spectrum 

Next, let us consider the time-dependent power spec- 
trum of the emitted light, which for sufficiently small 
passband width of the spectral apparatus can be 
given by [see, e.g., Vogel et al. (199^ )] 



S(t,us,T) 



1 f T 
dt 2 / dti 

o J o 



-iu>s{t2— ti) 



x(EH(r I t 2 )EW(r ) t 1 ))] 1 (3.38) 

where uj$ is the setting frequency of the spectral 
apparatus and T is the operating-time interval of 
the detector. In close analogy to the derivation of 
Eq. (|3~29|) , combination of Eqs. Q - Q, ( jD]| , 
Q , and Q leads to 



S(y,uj s ,T) 



T 



dti 



J(0J S -U] A )tl 



tl 



x / dt'Cu{t') / dwe^^-^X* 1 -''' 
Jo Jo 

2 



x ImG(r, Y A ,uj)d 



(3.39) 



Further calculation again requires knowledge of the 
Green tensor of the problem. 

Let us use Eq. ( |3.39|) to recover the free-space re- 
sult. Following the line that has led from Eq. (|3.2£ 
to Eq. ( 3.33j) , we find that 



S(r,us,T) 



( co A d sin 6 
\ 4ireoc 2 p 



e(T- P /c) 



-r /2+i(w s -uj' A )](T-p/c) _ x 



uj s -uj' a + To/2 



(3.40) 
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In particular for T^oo, we recognize the well-known 
Lorentzian: 



lim S(r,u>s,T) 



u> A dsm 6 

4:TTEoC 2 p 



(3.41) 



If retardation is ignored and the Markov approx- 
imation applies, Eq. ( |3.39|) can be simplified in a 



similar way as Eq. ( 3.29 ). In close analogy to the 
derivation of Eq. ( |3.34| ) we may write 



S(r,u s ,T) = \F(r,r A ,u A ) 



uis - uj' a + iT/2 
with ¥(r,r A ,u A ) from Eq. ( g3p . 



(3.42) 



3.3 Two- Atom Coupling 

We now turn to the problem of two atoms (denoted 
by A and B) coupled through a medium-assisted 
electromagnetic field in the case of single-quantum 
excitation. For simplicity, let us consider atoms with 
equal transition frequencies, so that 

K AA ,{tX) = K AA ,{t-t') (3.43) 

(A' = B,D) and 

K AB (t-t')=K BA (t-t'), (3.44) 

and assume that the isolated atoms undergo the 
same decay law, 

K AA (t,t') = K BB (t,t') = K(t-t'). (3.45) 

Introducing the new variables 



c±(t) 



-1/2 



C UA (t)±C UB (t) 



(3.46) 



it is not difficult to prove that the integro-differential 
equations (|3.10| ) decouple as follows: 



C±(t) = / dt' K ± (t - t') C±(0 
Jo 

+ 2- 1 ' 2 f dt' [K AD (t-t') 
J-At 

±K BD (t-t')]C UD (t'), (3.47) 

where 

K±(t - t') = K(t - t') ± K AB {t - t'). (3.48) 



Obviously, the C± are the expansion coefficients of 
the wave function with respect to the (atomic) basis 



2-V 2 (\U A )±\U B )), 



(3.49) 



and \L) (instead of the basis \U A ), \U B ), and \L)). 
Thus, they are the probability amplitudes of finding 
the total system in the states |+)|{0}) and |— )|{0}), 
respectively. In the further treatment of Eq. ( p.47[ ) 
one can again distinguish between the weak- and the 
strong-coupling regime. 



Weak Coupling In the weak-coupling regime, the 
Markov approximation applies, and in Eq. ( |3.47 ) 
C±(t') can be replaced with C±(t), with the time 
integrals being ("-functions. In particular, when the 
field is initially not excited, then the second term on 
the right-hand side of Eq. ( |3.47| ) vanishes and we are 
left with a homogeneous first-order differential equa- 
tion, whose solution is, in analogy to Eq. ( |3T2l| ) , 



C±(t) = e(- r± / 2+i5± )*C±(0), (3.50) 
where (r = T AA , 5 = 5 AA ) 

r± = r±I\4 B , (3.51) 



F AB 



5ab : 



heod 2 



V 



5± = 5±5 AB , (3.52) 
d A lmG(r A ,r B ,uj A )d B , (3.53) 
uj 2 d A lm G{r A ,r B ,uj)d B 



duj ■ 



o Jo 



u> — UJ A 



(3.54) 

Clearly, T± are the decay rates of the states |±), 
and the assumption ( |3.45| ) means that the two atoms 
are positioned in such a way that they have equal 
single-atom decay rates and line shifts. Note that 
the values of r + and T_ can substantially differ from 
each other, because of the interference term T AB (of 
positive or negative sign). 

Strong Coupling In the strong-coupling regime, 
the atoms are predominantly coupled (in a 
resonator-like equipment) to a sharp field reso- 
nance, whose mid-frequency approximately equals 
the atomic transition frequency. As a result, the 
atomic probability amplitudes in Eq. ( p.47[ ) must 
not necessarily be slowly varying compared with the 
kernel functions and the Markov approximation thus 
fails in general. Regarding the line shape of the field 
resonance as being a Lorentzian, one can of course 
approximate the kernels K(t — t'), K AB (t — t') [and 
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KAD(t — t') and KBD(t — t')] in a similar way as done 
in Eq. ( |3.24j) for a single atom. 

Equation ( 3.47] ) reveals that the motion of the 
states |±) defined by Eq. (|3.49| ) is governed by the 
kernel functions K±(t — tf), and it may happen that 
one of them becomes very small, because of destruc- 
tive interference [cf. Eq. ( |3.48| )1. In that case, either 
|+) or | — ) is weakly coupled to the field, and thus the 
strong-coupling regime cannot be realized for both 
of these states simultaneously. 

4 BULK MEDIUM 

The formalism outlined in Section |3| requires knowl- 
edge of the permittivity as a function of space and 
frequency. The spatial variation is typically deter- 
mined by some arrangement of macroscopic bodies, 
each of which being characterized by a permittivity 
that is a function of frequency only. The frequency 
response of a dielectric body depends on its atomic 
structure and can be measured with high precision. 
For theoretical studies it may be useful to have some 
analytical expression at hand. 

4.1 Drude-Lorentz Model 

In the Drude-Lorentz model, which is widely used 
in practice, the permittivity is given by 



Pa 



u). 



Ta 



where lOTa and j a are the medium oscillation fre- 
quencies and linewidths, respectively, and ujp a cor- 
respond to the coupling constants. It is worth noting 
that the Drude-Lorentz model covers both dielectric 
(wfa 7^ 0) and metallic {tOTa = 0) matter. An exam- 
ple of the permittivity for a (single-resonance) di- 
electric as a function of frequency is shown in Fig. |l[ 
From the permittivity, the refractive index can be 
obtained according to the relations 



n R (uj) + ini{u), 



(4.2) 



(4.3) 



The Drude-Lorentz model features band gaps be- 
tween the transverse frequencies UTa and the longi- 
tudinal frequencies u!L a = \/uj, 



— / c ~ 


i 

- 







0.9 



(D/(D T 



1.1 



Figure 1: Real and imaginary parts of the permit- 
tivity of a single-resonance Drude-Lorentz-type di- 
electric for cop = 0.5 ujt and 7=10~ 4 wt- The band 
gap covers the interval from ujt to ujl — 1 -12 lot- 

medium resonance, we typically observe that 

eiiu) < \e R {u)\. (4.4) 

For ijj < u>Ta (outside a band gap) we have 

e r {u) > 1, (4.5) 

n R {uj) ~ y/e R {u) > m{uj) ~ - - — , 

2 v / e jR ,(w) 

(4.6) 



(4.1) and for ujt q < w < u>ia (inside a band gap) 



e R (u) < 0, 



n R {u)) 



e/(w) 



2y^H| 



< ni(u) ~ yJ\e R {uo) \ . 



When (inside a band gap) 

e R {u) < -1 
is valid, which, in view of Eq. ([4.1|), leads to 



u < 



U Ta + 2 U Pa ' 



(4.7) 

(4.8) 
(4.9) 

(4.10) 



then the Drude-Lorentz model also incorporates 



surface-guided waves [see, e.g., Raether (1988 ); Ho 
et al. (20011) 1, which are observed in the presence 



Ta 



+ UJ 



Pa- 



Far from a 



of an interface. These waves are bound to the in- 
terface, with the amplitudes being damped into ei- 
ther of the neighboring media. Typical examples are 
surface phonon polaritons for dielectrics and surface 
plasmon polaritons for metals. Note that in any case 
ei(u) > is valid. 
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4.2 Local-Field Correction 



Fro' in simple arguments based un bfae change uf bfae 

mode density, it was suggested that the sponta- 
neous emission rate of an atom inside a nonabsorbing 
medium should be modified according to T = nTo, 
where n is the (real) refractive index of the medium 
and Tq is given by Eq. Q3.18 ) [ Dexter (1956| ); |di Bar- 
tolo (1968| ); [Yariv (1975| ); [Nienhuis et al. (1976Q ]. 
In this formula, it is assumed that the local field 
the atom interacts with is the medium-assisted elec- 
tromagnetic field obtained by averaging over a re- 
gion which contains a great number of medium con- 
stituents. In reality, the atom is located in a small 
region of free space, and the field therein differs from 
the averaged field. This effect is usually taken into 
consideration by introduction of a local-field correc- 
tion factor £, thus 



[[Rikkcn et al. (1995| ); |de Vries et al. (1998Q ; (Schuur- 
mans et al. (1998|) ] . 



<r . 



(4.11) 



In the (Clausius-Mosotti) virtual-cavity model it is 
given by [[Knoester et al. (1989| ); [Milonni (19951 )1 



n 2 + 2 



(4.12) 



and in the (Onsager) real-cavity model by | Glauber 



et al. (199 ID 



3n 



2 \ 2 



2n 2 + l/ 



(4.13) 



For absorbing media it was suggested that the in- 
dex of refraction should be replaced by its real part 
and the square of the correction facto r in Eqs. ([4.1 2j) 
and ( 4.13 ) by the absolute square [ Barnett et al. 
(1992| , |1996|) ; [luzeliunas (19971 )1- Later it was found, 



on the basis of the quantization scheme outlined in 
Section |2[ that a proper inclusion of the (quantum) 
noise polarization leads to a more complicated form 
of the local-field correction [[Scheel et al. (1999a| ,[b|)1, 



which (for the virtual-cavity model) was confirmed 



by an alternative, microscopic approach | Fleisch- 
hauer (19991 )]. 



In contrast to the virtual cavity model where the 
modification of the field outside the cavity is disre- 
garded, in the real-cavity model the mutual mod- 
ification of the field outside and inside the cavity 
are taken into account in a consistent way. Exper- 
iments suggested that the real-cavity model is suit- 
able for describing the decay of substitutional guest 
atoms different from the constituents of the medium 



Let us consider an excited atom placed at the 
center of an empty, spherical cavity (embedded in 
an otherwise homogeneous medium). According to 
Eq. ( |3.22j ) and the Green tensor given in Appendix 
[A.2| , the decay rate can be given in the form of 
[|5cheel et al. (1999b|)] 



r = r [i + ReCf (u A )] 



(4.14) 



where the generalized reflection coefficient C± \u) 
reads \k = k{ui) = Rlo/c] 

Cf (u) = e fk {i + fc[n(w)+l] - ik 2 n{u) 

- k 3 n 2 (lo) I [n(u) + 1] } { sin k - k [cos k 
+ in(uj) sin k]+ik 2 n(uj) cos k — k 3 [cos k 

— in{uj) sin k}n 2 (cj)/[n 2 (uj)-l]} 1 (4.15) 

As long as the surrounding medium can be treated 
as a continuum, Eq. ( 4.14 ) [together with Eq. ( [4.15 )1 
is exact. It is valid for arbitrary cavity radius and 
arbitrary complex refractive index, without restric- 
tion to transition frequencies far from medium reso- 
nances. 

When the cavity radius is much smaller than the 
wavelength of the atomic transition, i.e., Rua/c = 
k(uA) *C 1, then the real-cavity model of local-field 
correction is realized. In this case, {uja) can be 
expanded in powers of k{uJA) to obtain [ [Schccl et al. 
(1999bD ; [Tomas (200l| )1 



r = r n 



3e 



2e+l 



+ 



£l 



< \ 3 + 28|e| 2 + ](i: ;r - 1 / c 



loaR 



5|2e + l| 2 \u;aR 



|2e+l| 2 
+ 0{u A R/c), 



2ni\e\ + niER+n R Ei) 



(4.16) 



where the dependence of the permittivity e and the 
refractive index n on loa has been suppressed. For 
ei(uja) = 0, i.e, when material absorption is fully dis- 
regarded, Eq. ( 4.16[) r eproduces exactly local-field 
correction factor ( 4.13| ). The second term in the 
curly brackets essentially results from absorption. It 
is seen that material absorption gives rise to a strong 
dependence of the decay rate on the cavity radius. In 
particular, the leading term proportional to R~ 3 can 
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half-spaces such that 



e\r,u 



1 



if z < 
if z > 



(5.1) 



For z > 0, the reflection part of the Green tensor 
reads fMaradudin et al. (1975| ); [Mills et al. (1975|) ; 
Tomas (1995| ); |Ho e* a/. (19981) ] 



Gx X (z, z, u) 



irk 2 j 

dfc„^^r s (fc||), 



dfc|| fc||/3e^ 2 r p (/fc||) 



(0 A /(0 7 



+ ■ 



Figure 2: Spontaneous decay rate [Eq. ([4,14|) 1 
of an atom embedded in a dielectric medium as 
a function of the atomic transition frequency near 
a medium resonance for a single-resonance Drude- 
Lorentz-type dielectric [ujp = 0.5ujt] "//uJt = 10~ 2 
(solid line), 10 -3 (dashed line), and 10~ 4 (dotted 
line)]. The (real-) cavity radius is R = 0.02 Xt- 



be regarded as corresponding to nonradiative energy 



8vrjo 

Gyy(z, Z, UJ 



1^ 



Gx X (z, z, uj), 



4nk 2 



(5.2) 
(5.3) 
(5.4) 



tra nsfer from the atom to the medium. Examples 



[k = uj/c, (3 = (k 2 — k 2 ) 1 / 2 ], where r p (k\\) and r s (k\\) 
are respectively the familiar Fresnel reflection coef- 
ficients for p- (TM) and s- (TE) polarized waves. 
Substitution of these expressions (uj = uja) into 
Eqs. ( |3.22| ) and ( p. 23 ) yields the decay rate and the 
line shift of an atom at a distance z from the sur- 
face [ Agarwal (1975 ); Agarwal et al. (1977 ); (School 
et al. (1999c| )], which are in agreement with classical 



of the dependence of rate of spontaneous decay on 
the atomic transition frequency are plotted in Fig. |2| 
for a Drude-Lorentz-type dielectric medium. It is 
seen that in the band-gap region (where for a non- 
absorbing medium spontaneous emission would be 
inhibited) the decay rate can drastically increase, 
because of the non-radiative decay channel associ- 
ated with absorption. Note that the strongest en- 
hancement of spontaneous decay is observed at uja 

uj\-\-\uj 2 p , which [for small values of ei(uja)] cor- 



results [see, e.g., |Chance et al. (1978| ) and references 
therein] . 

When the distance of the atom from the surface 
is small compared to the wavelength, b<l, then 
the integrals in Eqs. (|5.2|) - (|5.4|) can be evaluated 
asymptotically to give [{Bcheel et al. (1999^ )] 

n 2 (uj) 



G? z (z,z,uj) 



1 



1 



responds to 2e(uja) + 1 — 0. 



5 PLANAR SURFACE 

Let us turn to the problem of spontaneous decay of 
an excited atom located near the surface of a half- 
space medium. For real permittivity, configurations 
of that type have been studied extensively in con- 
nection with Casimir and van der Waals forces [see, 

Fichet et al. (1995| ) 



+ 



+ 



1 



n[uj 



16irk 2 z 3 n 2 (uj) + 1 
)-l] 2 



8ttz n(uj)[n(uj) + 1] 
ik [n(oj) - l][2n(w) 



1] 



12?r n(u)[n(u) + 1] 



Gxx ( z i z i w ) — \G^ z (z, z, uj) - 



ik n(uj) 



Q(kz), (5.5) 



1 n 2 (uj)-l 



16nz n 2 {uj) + l 



- + 0(kz), 
3ir n(uj) + 1 

Gy V (z,z,uj) = Gxx(z,z,uj). 



e.g., [Meschede et al. (1990j) 



and references therein] and with regard to scanning 



Inserting Eqs. (|5j| - (|J) into Eq. (|3.2 

3 



(5.6) 

(5.7) 
yields 



3r, 



near-field optical microscopy [see, e.g., Henkel et al. 



1 -I — £ 
+ d 2 



uj a z 



(19981) 1 . 



To be more specific, let us consider two infinite 



£j(^a) 

k(^) + i| 2 



UJ A Z 



(5.8) 
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The leading (~ z -3 ) term is the same as in the micro- 
scopic approach by Yeung et al. (1996| ). This term, 
which is proportional to ei(uja), is closely related 
to nonradiative decay, i.e., energy transfer from the 
atom to the medium. Obviously, a change of ei{uja) 
mostly affects the near-surface behavior of the decay 
rate. Note that the distance of the atom from the 
surface must not be smaller than interatomic dis- 
tances in the medium (otherwise an interface can- 
not be defined). Examples of the spontaneous de- 
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Figure 3: The rate of spontaneous decay of 
an excited atom near a planar dielectric half- 
space is shown as a function of the transition fre- 
quency for an x-oriented transition dipole moment 
and a single-resonance Drude-Lorentz-type dielec- 
tric [ujp = 0.5 ujt; 7/k>r = lCr 2 (solid line), 1CT 3 
(dashed line), and 10~ 4 (dotted line); z = 0.05 At]- 
The inset illustrates the dependence on the distance 
of the decay rate [j/ljt = 10~ 3 ; ooa/^t = 0.9 (solid 
line), 1.06 (dashed line)]. 

cay rate are shown in Fig. || for dielectric matter of 
Drude-Lorentz type. Note the strong absorption- 
assisted enhancement of spontaneous decay that is 

observed inside the band gap at uja — yj ui^ + ^ojp. 
It corresponds [for small values of £i(u)a)] to the 
condition s(uja) + 1 — 0, which marks the position of 
the highest density of the surface-guided waves [cf. 
Eq. O]. 



Similarly, from Eq. (3.23) together with Eqs. (|5^ 



- (5.7), the line shift due to the presence of the 
macroscopic body reads 



6uja 



3To 
32 



1 + 



4 
i 



LO A Z 



H"a) + 1\ 2 



+ o[ 



\UJ A Z 



(5.9) 



In contrast to the decay rate, here the leading 
(~ z~ 3 ) term even appears when absorption is disre- 
garded [ei(u A ) = 0]. 

6 SPHERICAL MICRO- 
RESONATOR 



In Section 4.2, an atom in a microsphere whose ra- 
dius is much smaller than the wavelength of the 
atomic transition was considered. If the radius is 
not small compared with the wavelength, the cav- 
ity can act as a resonator. It is well known that 
the spontaneous decay of an excited atom can be 
strongly modified when it is placed in a microres- 
onato r [[Hinds (1991|); |Haroche (199% [Meschede 
" 19921 ); |Meystre (1992|) ; |Berman (Ed., 1994j) ; |Kim- 
There are typically two qualitatively 
the weak-coupling regime and 



ble (1998| ; 
different regimes: 
the strong-coupling regime. In the weak-coupling 
regime the Markov approximation applies and a 
monotonous exponential decay is observed, the de- 
cay rate being enhanced or reduced compared to the 
free-space value depending on whether the atomic 
transition frequency fits a cavity resonance or not. 
The strong-coupling regime, in contrast, is charac- 
terized by reversible Rabi oscillations where the en- 
ergy of the initially excited atom is periodically ex- 
changed between the atom and the radiation field. 
This usually requires that the emission is in reso- 
nance with a high-quality cavity mode. 




Figure 4: Scheme of the spherical microresonator. 

Let us consider an excited atom placed at the cen- 
ter of a spherical three-layer structure (Fig. ||). The 
outer layer (r > R\) and the inner layer (0 < r < R2) 
are vacuum, whereas the middle layer {R2 < r < 
Ri), which plays the role of the resonator wall, is 
matter. In particular for a Drude-Lorentz-type di- 
electric, the wall would be perfectly reflecting in the 
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band-gap zone, provided that absorption could be 
disregarded. Restricting our attention to a true res- 
onator, we may assume that the condition B,2Uja/c^> 
1 is satisfied. 

6.1 Weak Coupling 



From Eq. (1.3) together with the Green tensor for a 
spherical three-layer structure as given in Appendix 
A. 2, the decay rate becomes JHo et al. (200C| )] 



r Re 



n(u A ) — i ta.n{uj aR2 / c) 



1 — in(u)A) tan(a>A-R2/c) 
r n R (uj A )[l + taxi 2 (lu a R 2 /c)} 

x | [1 + tii((jJa) tan(uJAR2/c)] 2 
+ u r(. uj a) tan 2 (aj A R 2 /c)} . 



(6.1) 



Note that in Eq. (p.l|) it is assumed that 
exp[— ini(coA)(Ri — R2)^a/c] <1 (thick cavity 
wall) . 
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inhibition (T/Tq < 1). The frequencies at which the 
maxima of enhancement are observed correspond to 
the resonance frequencies of the cavity. Within the 
band gap the heights and widths of the frequency 
intervals in which spontaneous decay is feasible are 
essentially determined by material absorption. Out- 
side the band-gap zone the change of the decay rate 
is less pronounced, because of the relatively large 
input-output coupling, the (small) material absorp- 
tion being of secondary importance. 




C0 A /C07 



Figure 6: The amount of the (outside the resonator) 
available radiation energy after spontaneous decay 
of an excited atom at the center of a spherical mi- 
croresonator is shown as a function of the transition 
frequency for a single-resonance Drude-Lorentz-type 
dielectric wall [j/ujt = W~ 2 (solid line), 2 x 1CT 2 
(dashed line), 5 x 10 -2 (dotted line); the other pa- 
rameters are the same as in Fig. EJ. [After Ho et al 



Hoot]).] 



Figure 5: The rate of spontaneous decay of an 
excited atom at the center of a spherical microres- 
onator is shown function of the transition 

frequency for a single-resonance Drude-Lorentz- 
type dielectric wall (i?2 = 30Ar; R\ — R2 = ^t] 
wp = 0.5u;t; 7 = 10 _2 u>r)- The curves in the in- 
set correspond to ^/ujt = 10~ 2 (solid line), 2 x 1CT 2 
(dashed line), and 5 x 1CT 2 (dotted line). [After |Ho| 
\et al. (2000Q .] 

The dependence on the transition frequency of the 
decay rate is illustrated in Fig. |5[ It is seen that the 
decay rate very sensitively depends on the transi- 
tion frequency. Narrow-band enhancement of spon- 
taneous decay (T/Tq > 1) alternates with broadband 



The widths of the resonance lines are responsi- 
ble for the damping of the corresponding intracavity 
field (mode). There are two damping mechanisms: 
photon leakage to the outside of the cavity and pho- 
ton absorption by the cavity-wall material. The 
first mechanism is the dominant one outside bands 
in regions where normal dispersion (dnn/duj > 0) is 
observed, while the latter dominates inside band 
gaps where anomalous dispersion (dn^/dw<0) is 
observed. To illustrate this in more detail, let us 
consider the total amount of radiation energy ob- 
served outside the cavity and compare it with the 
energy Wq = Huja emitted by an atom in free space. 
Application of Eq. ( 3.37 ) [together with Eqs. ( |3.34 ) 
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and ( pT36D ] yields [|Ho et al. (20"00D ] 



W 

Wo 



l + ReC/V A ) ' 



(6.2) 



with ^^(wa 



and Cf( 



ojx) being given according to 
Eqs. ( A. 34 ) and ( A.35|) respectively. Examples of 
the dependence of W/Wq on the atomic transition 
frequency are plotted in Fig. ^. It is seen that inside 
the band gap most of the energy emitted by the atom 
is absorbed by the cavity wall in the course of time, 
while outside the band gap the absorption is (for the 
chosen values of 7) much less significant. Note that 
with increasing value of 7 the band gap is smoothed 
a little bit, and thus the fraction of light that escapes 
from the cavity can increase. 

6.2 Strong Coupling 

When the coupling between the atom and a cavity 
resonance (mid- frequency loc, line width Aloc) is so 
strong that the (weak-coupling) decay rate becomes 
comparable to the cavity line width, Tc <; Aujc, the 
Markov approximation is no lon.qpr adequate Tn this 



cas e, the integral equation ( 3rjrS ) or approximate 



equations of the type ( |3.25|) and ( |3.26 ) should be 
used in order to describe the temporal evolution of 
the (upper) atomic state. For the configuration un- 
der investigation, the cavity line width can be given 
by 

cT 



Awe 



(6.3) 



R2T C 

Inside a band gap, Tc is essentially determined 
by material absorption. In particular, the single- 
resonance Drude-Lorentz model reveals that 



C 



r, 



nj(uj c ) + 1 
Ur{u)c) 



~ r 



uJ. 



r) 



7^c 



(6.4) 



(7 Cwt,wp,Wp/wt)- Below the band gap radiative 
losses dominate and Tc reads 



c 



To n R (u;c) 











(jjj, 


-uj 2 c 



(6.5) 



(nR>n/). 

Typical examples of the time evolution of the 
upper-state occupation probability are shown in 
Fig. The curves are the exact (numerical) so- 
lutions of the integral equation ( |3.13j ) [together with 




Figure 7: The temporal evolution of the occupation 
probability of the upper state of an initially excited 
atom at the center of a spherical microresonator is 
shown for a single-resonance Drude-Lorentz-type di- 
electric wall [R2 = 30 At; R\ — R2 = At; ^>p = 0.5 ojt\ 
u A = 1.046448 u T ; T X T /(2c) = 10~ 6 ; 7/w r = 10- 4 
(solid line), 5 x 10 -4 (dashed line), 10~ 3 (dotted 
line)]. For comparison, the exponential decay in free 
space is shown (dashed-dotted line) . [After [Ho et al . 

B )-1 



the kernel function ( 3.17] )] for a (single-resonance) 
dielectric wall of Drude-Lorentz type. The figure 
shows that with increasing value of the intrinsic ab- 
sorption constant 7 of the wall material the Rabi 
oscillations become less pronounced. Clearly, larger 
values of 7 mean enlarged absorption probability of 
the emitted photon by the cavity wall and thus re- 
duced probability of atom-field energy interchange. 



7 MICROSPHERE 

Light propagating in a dielectric sphere can be 
trapped by repeated total internal reflections. When 
the round-trip optical path fits integer numbers of 
the wavelength, whispering gallery (WG) waves are 
formed, which combine extreme photonic confine- 



ment with very high quality factors [Collot et al. 



19931 ); |Chang et al. (1996|) ; |Gorodetsky et al. (1996|) 



Vernooy et al. (1998| ); [Uetake et al. (1999j )1 - prop- 



erties that are crucial for cavity QED experiments 
[ [Lin et al. (1992 ); Barnes et al. (1996| ); Lermer et al. 



T99l|); |Vernooy et al (1998| ); [Fujiwara et al. (1999| ) 
Yukawa et al. (1999| )1 and certain optoelectronical 
applications [Chang et al. (1996)]. WG waves are 



|commonly classified by means of three numbers |Col- 



1 



lot et al. (1993); |Chang et al. (19961 )]: the angular- 
momentum number I, the azimuthal number m, and 
the number i of radial maxima of the field inside the 
sphere. In the case of a uniform sphere, the WG 
waves are (2Z + l)-fold degenerate, i.e., the 21 + 1 
azimuthal resonances belong to the same frequency 

A dielectric microsphere whose permittivity is of 
Drude-Lorentz type does not only give rise to WG 
waves, but can also feature surface-guided (SG) 
waves inside band-gap regions.^ In contrast to WG 
waves, each angular-momentum number / is associ- 
ated with only one SG wave. 

If an excited atom is situated near a dielectric mi- 
crosphere, spontaneous decay sensitively depends on 
whether or not the transition is tuned to a WG or an 
SG resonance. Moreover, whereas WG waves typi- 
cally suffer from material absorption, the effect of 
material absorption on SG waves is weak in general. 

7.1 Decay Rate 



100 



Applying Eq. ( |3.22j ) together with the Green tensor 
of a microsphere (Appendix A. 2 ), the spontaneous 
decay rate for a (with respect to the sphere) radially 
oriented dipole moment can be given by 



r ± = r i + lE 



i=i 



/(/ + 1)(2/ + 1) 



and for a tangential dipole it reads 



(7.1) 



i=i 



(21 + 1) 



xRefeVA) 



hf\k A r A ) 



[kATAh^ {k A rA)]' 



kATA 



(7.2) 



(kA = uja/c), where B t , (wyi) is defined according to 
Eq. (|A.30|) , and the prime indicates the derivative 
with respect to kA^A [ Ho et al. (20Q1|) ] .P| Note that a 



For the dependence on frequency of the quality factors of 



Ho et al. (2001) 



WG and SG waves, see 

3 Equations (7.1) and (7.2) can be regarded as being the 



natural extension of the (clas sical ) theory for nonab sorbing 



matter as given by Shew (1987 



and 



Klimov et al. (1996). Note 



that when the formulas for nonabsorbing matter are given in 
terms of the Green tensor, without an explicit decomposition 
in real and imaginary parts, then for the real permittivity the 
complex one can be substituted. 




Figure 8: The rate of spontaneous decay of an 
excited atom near a microsphere is shown ELS cl 
function of the transition frequency for a radially 
oriented transition dipole moment and a single- 
resonance Drude-Lorentz-type dielectric [R = 2Xt', 
ujp = 0.5 ujt', 7/wt = 10 -4 ; (a) Ar = 0.02Ar; inset: 
7/Vr=l(T 4 (solid line), l(T 5 (dashed l ine), 10~ 6 
(dotted line); (b) Ar = 0.1A T ]. [After |Ho et al 
T200l )-1 



radially oriented transition dipole moment only cou- 
ples to TM waves, whereas a tangentially oriented 
dipole moment couples to both TM and TE waves. 

It is worth noting that when the atom is very close 
to the microsphere, then the decay rates Eqs. (|7. 1[) 
and (|7.2| ) reduce to exactly the same form as in 
Eq. ( |5.8[) for a planar interface, with z being now 
the distance between the atom and the surface of the 
microsphere. Obviously, nonradiative decay, which 
dominates in this case, does not respond sensitively 
to the actual radiation-field structure. 

The dependence on the transition frequency of 
the decay rate, as it is typically observed for not 
too small (large) values of the atom-surface distance 
(material absorption), is illustrated in Fig. || for a 
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radially oriented transition dipole moment. Since 
it mimics the single-quantum excitation spectrum 
of the sphere-assisted (TM) radiation field, the fig- 
ure reveals that both the WG and SG field reso- 
nances can strongly enhance the spontaneous decay. 
Material absorption broadens the resonance lines at 
the expense of the heights, and the enhancement 
is accordingly reduced [see the inset in Fig. |8](a)]. 
Clearly, the sphere-assisted enhancement of sponta- 
neous decay decreases with increasing distance be- 
tween the atom and the sphere [compare Figs. |8|(a) 
and (b)]. 

Figure || also reveals that SG waves can give rise 
to a much stronger enhancement of the spontaneous 
decay than WG waves. In particular, with increasing 
angular-momentum number the SG field resonance 
lines strongly overlap and huge enhancement [e.g., 
of the order of magnitude of 10 4 for the parameters 
chosen in Fig. 0(a)] can be observed for transition 
frequencies inside a band gap. When the distance 
between the atom and the sphere increases, then the 
atom rapidly decouples from that part of the field. 
Thus, the huge enhancement of spontaneous decay 
rapidly reduces and the interval in which inhibition 
of spontaneous decay is typically observed, extends 
accordingly [see Fig. ||(b)]. 

7.2 Frequency Shift 

The sphere-assisted frequency shift calculated from 
Eq. ( |3.23|) together with the Green tensor given in 
Appendix [A. 2 reads 
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for a radially oriented transition dipole moment, and 
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(7.4) 



for a tangentially oscillating dipole.0 Note that the 
small quantum corrections that arise from the sec- 



4 Again, Eqs. ( |7.3[ ) and (7.4) could be obtained from the 
classical theory for nonabsorbing matter [Klimov et al. (1996)]; 
see footnote <A on page 16. 
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Figure 9: The frequency shift in spontaneous 
decay of an excited atom near a microsphere 
is shown as a function of the transition fre- 
quency for a radially oriented transition dipole mo- 
ment and a single-resonance Drude-Lorentz-type di- 
electric [R = 2Xt; u>p = 0.5u>t', 7/wt = 10 -4 ; (a) 
Ar = 0.02A T ; (b)Ar = 0.1A T ; inset: 7/cj t =10" 4 
(solid line), 10~ 5 (das hed line), 10 -6 (dotted line)]. 
[After |Ho et al. (200l] ).l 



ond term in Eq. (3.23) have been omitted. For very 
small distance between the atom and the sphere, 
Eqs. (|7.3j ) and ([M|) acquire the same form as for 
a planar interface, Eq. ( |5.9| ). 

In Fig. ^, an example of the dependence on the 
transition frequency of the frequency shift for a ra- 
dially oriented dipole is shown. It is seen that the 
field resonances can give rise to noticeable frequency 
shifts in the very vicinities of the corresponding 
resonance frequencies. Transition frequencies that 
are lower (higher) than a resonance frequency are 
shifted to lower (higher) frequencies. In close anal- 
ogy to the behavior of the decay rate, the frequency 
shift is more pronounced for SG resonances than for 
WG resonances and can be huge for large angular- 
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momentum numbers when the lines of the SG field 
resonances strongly overlap. 

The behavior of the frequency shift as shown in 
Fig. ^(b) can already be seen in the single-resonance 
approximation [ |Ching et al. [1987] )] . Let the atomic 
transition frequency uj a be close to a resonance fre- 
quency u>c of the microsphere and assume that, in a 
first approximation, the effect from the other reso- 
nances may be ignored. For a Lorentzian resonance 
line of width Auc, from Eq. Ql.4| ) it then follows 
that 



V 

47T 



U> — u>a (uJ—UJc) 2 +&-u'c 



UJ A - UC 



2 (lo a - lo c ) 2 + Acj£ 



2 ' 



(7.5) 



where Tc (which corresponds to the height of the 
line) is the decay rate for loa = ^>c- In particular, 
Eq. ( |7.5|) indicates that the frequency shift peaks at 
half maximum on both sides of the resonance line. 
With increasing material absorption, the linewidth 
Au>c increases while Tc decreases, and thus the ab- 
solute value of the frequency shift is reduced, the 
distance between the maximum and the minimum 
being somewhat increased. With decreasing distance 
between the atom and the microsphere near-field ef- 
fects become important and Eq. (|7.5| ) fails, as it can 
be seen from a comparison of Figs. |9|(a) and (b). 

7.3 Emitted-Light Intensity 
7.3.1 Spatial distribution 

Substitution into Eq. ( |3.35|) of the expression for the 



Green tensor (Appendix A.2 ) yields (9a = (Pa = 0, 

ta < r ) 
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for a radially oriented transition dipole moment, and 
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for a tangentially oriented dipole in the xz-plane. 
Here the abbreviating notations 



B 



N 



1 



kATA 



[kArAjl(kArA)}' 



- Bf(u A ) [k A rAhf ] (k A r A )}'}, (7.8) 
Bf = ji(k A r A ) + Bf*(u>A)h?\kArA), (7.9) 

P^cosO) = 1(1 + l)P,(cos#) - cos9P((cos9) (7.10) 

have been introduced. |F- L tl^ (r, r A , lva)\ 2 deter- 
mines, according to Eq. ( |3.34D , the spatial distribu- 
tion of the light emitted by a radially (tangentially) 
oriented dipole. 

Let us restrict our attention to a radially ori- 
ented transition dipole moment. Examples of 
|F J -(r, ta, ua)\ 2 are plotted in Fig. [H]. In this case, 
the far field is essentially determined by Fjj~, as an 
inspection of Eq. fl7.6j ) reveals. When the atomic 
transition frequency coincides with the frequency of 
a WG wave of angular-momentum number I far from 
the band gap [Fig. |i~0|(a )], then the corresponding 
Z-term in the series ( |7.6|) obviously yields the lead- 
ing contribution to the emitted radiation, whose an- 
gular distribution is significantly determined by the 
term ~sin#i^'(cos#). Thus, the emission pattern 
has I lobes in, say, the yz-plane, i.e., I cone-shaped 
peaks around the z-axis, because of symmetry rea- 
sons. The lobes near 9 = and 9 = tt are the most 
dominant ones in general, because of 

-sin0P/(cos0) ~ (sin 9)' 1 /2 + 0(r 1 ) (7.11) 

(0 < 9 < it). Note that the superposition of the lead- 
ing term with the remaining terms in series ( |7.6| ) 
gives rise to some asymmetry with respect to the 
plane 9 = n/2. 

When the atomic transition frequency approaches 
(from below) a band gap (but is still outside it), a 
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Figure 10: Polar diagrams of the far-field emission pattern |F (r, r^, uja)\ 2 / {k\d/ Atteq) 2 of a radially 
oscillating dipole near a single-resonance Drude-Lorentz-type dielectric microsphere [R = 2\t; ujp = 0.5lvt] 
7/ Vr = l (T 4 ; Ar = 0.02A T ; r = 20A T ; uj a /oj T = 0.94042 (a), 0.999 (b), 1.02811 (c), 1.06 (d). [After [Ho~^ 
al. [2001 ). 1 



strikingly different behavior is observed [Fig. 10(b)]. 
The emission pattern changes to a two-lobe structure 
similar to that observed in free space, but bent away 
from the microsphere surface, the emission intensity 
being very small. Since near a band gap absorption 
losses dominate, a photon that is resonantly emitted 
is almost certainly absorbed and does not contribute 
to the far field in general. If the photon is emitted in 
a lower-order WG wave where radiative losses domi- 
nate, it has a bigger chance to escape. The superpo- 
sition of all these weak (off-resonant) contributions 
just form the two-lobe emission pattern observed, 
as it can also be seen from careful inspection of the 
series ( |7.6j ). 

When the atomic transition frequency is inside 
a band gap and coincides with the frequency of a 
SG wave of low order such that the radiative losses 
dominate, then the emission pattern resembles that 
observed for resonant interaction with a low-order 
WG wave [compare Figs. 0(a) and (c)]. With in- 
creasing transition frequency the absorption losses 
become substantial and eventually change the 



emission pattern in a quite similar way as do below 
the band gap [compare Figs. [U](b) and (d)]. Obvi- 
ously, the respective explanations are similar in the 
two cases. 



7.3.2 Radiative versus nonradiative decay 

Since the imaginary part of both the vacuum Green 
tensor G v and the scattering term G R is transverse, 
the decay rate ( |1.3[ ) results from the coupling of the 
atom to the transverse part of the electromagnetic 
field. Nevertheless, the decay of the excited atomic 
state must not necessarily be accompanied by the 
emission of a real photon, but instead a matter quan- 
tum can be created, because of material absorption. 
To compare the two decay channels, let us consider, 
according to Eq. (|3.37 ), the fraction W/Wq of the 
atomic (transition) energy that is irradiated by an 
atom with a radially oriented transition dipole mo- 
ment. Using Eqs. ( 3.34 ) and ([Fl]), one derives [ po| 
et al. (20011 )1 
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Figure 11: The fraction of emitted radiation energy 
in spontaneous decay of an excited atom placed near 
a microsphere is shown as a function of the transition 
frequency for a radially oriented transition dipole 
moment and a single-resonance Drude-Lorentz- 
type dielectric [R = 2Xt; ujp = 0.5u>t', 7/wt = 10~ 4 ; 
Ar = 0.02 A T ; 7/w t = 10~ 4 (a), 1(T 6 (b). [After [Ho 
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Recall that W/Wq = 1 implies fully radiative decay, 
while W/Wq = implies fully nonradiative one. 

The dependence of the ratio W/Wq on the atomic 
transition frequency is illustrated in Fig. 11. The 



minima at the WG field resonance frequencies indi- 
cate that the nonradiative decay is enhanced relative 
to the radiative one. Obviously, photons at these 
frequencies are captured inside the microsphere for 
some time, and hence the probability of photon ab- 
sorption is increased. For transition frequencies in- 
side a band gap, two regions can be distinguished. 
In the low-frequency region, where low-order SG 



waves are typically excited, radiative decay domi- 
nates. Here, the light penetration depth into the 
sphere is small and the probability of a photon be- 
ing absorbed is small as well. With increasing atomic 
transition frequency the penetration depth increases 
and the chance of a photon to escape drastically di- 
minishes. As a result, nonradiative decay dominates. 
Clearly, the strength of the effect decreases with 
decreasing material absorption [compare Fig. [TT|(a) 
with (b)]. 

From the figure two well pronounced minima of 
the totally emitted light energy, i.e., noticeable max- 
ima of the energy transfer to the matter, are seen for 
transition frequencies inside the band gap. The first 
minimum results from the overlapping high-order 
SG waves that mainly underly absorption losses. 
The second one is observed at the longitudinal reso- 
nance frequency of the medium. It can be attributed 
to the atomic near-field interaction with the longitu- 
dinal component of the medium-assisted electromag- 
netic field, the strength of the longitudinal field res- 
onance being proportional to ej. Hence, the dip at 
the longitudinal frequency of the emitted radiation 
energy reduces with decreasing material absorption 
and may disappear when the atom is moved suffi- 
ciently away from the surface. 

7.3.3 Temporal evolution 

Throughout this section we have restricted our at- 
tention to the weak-coupling regime where the ex- 
cited atomic state decays exponentially, Eq. ( p. 21 ). 
When retardation is disregarded, then the intensity 
of the emitted light (at some chosen space point) 



simply decreases exponentially, Eq. (|3.34| ) . To study 
the effect of retardation, the frequency integral in 
the exact equation ( |3.29| ) must be performed numer- 
ically in general. 

Typical examples of the temporal evolution of the 



far-field intensity are shown in Fig. 12 for a radi- 
ally oriented transition dipole moment in the case 
when the atomic transition frequency coincides with 
the frequency of a WG wave. Whereas the long- 
time behavior of the intensity of the emitted light 
is, with little error, exponential, the short-time be- 
havior (on a time scale given by the atomic decay 
time) sensitively depends on the quality factor [Q 
~ 10 3 in Fig. 0(a), Q ~ 10 4 in Fig. 0(b)]. The 
observed delay between the upper-state atomic pop- 
ulation and the intensity of the emitted light can be 
quite large for a high-Q microsphere, because the 
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Figure 12: Exact (solid lines) and approximate 
(dashed lines) temporal evolution of the far-field 
intensity I(r, t)/(k A d/A7reo) 2 at a fixed point of 
observation of a radially oscillating dipole near 
a single-resonance Drude-Lorentz-type dielectric 
microsphere [R = 2Xt; tup = 0.5 u>t; 7/wt = 10 -4 ; 
Ar = 0.02A T ; r = 20A T ; 9 = 3; T /ojt = W 7 ; 
uj a = 0.91779 oj t (a 
20011 ). 1 



Figure 13: The decay rate of an atom near a 
metallic microsphere is shown ELS cl function of the 
transition frequency for a radially oriented transi- 
tion dipole moment and a single-resonance metal 
of Drude type [R = 5X P ; -//uj p = 0.005; Ar = 0.1A P . 
Inset: polar diagram of the far-field emission pat- 
tern |F J -(r, ta, ll>a)\ 2 /(k^fi/Aneo) 2 for r = 50Ap and 



al 



0.94042 uj t (b). [After [Ho ej u A /up = 0-5026. [After |Ho e* qZ. (200T| ).] 



time that a photon spends in the sphere increases 



with the (Rvalue. Further, m the short-time domain 
some kink-like fine structure is observed, which ob- 
viously reflects the different arrival times associated 
with multiple reflections. 

7.4 Metallic Microsphere 

The permittivity of a metal (on the basis of the 
Drude model) can be obtained by setting in Eq. Q4.1| ) 
the lowest resonance frequency loto equal to zero. 
Thus, the results derived for the band gap of a dielec- 
tric microsphere also apply, for appropriately chosen 
values of the corresponding ujp a = cop and j a = 7, to 
a metallic sphere. In particular, the results obtained 
in the nonretardation limit (c^oo) and for small 
sphere sizes (i?<A P ) URuppin (1982|) ; |Agarwal ej 
al . (1983| )1 are recovered. 



The dependence of the decay rate on the transi- 
tion frequency of an excited atom near a metallic 



microsphere is illustrated in Fig. 13 for R> Xp. The 
inset shows the emission pattern for the case when 
the atomic transition frequency coincides with the 
frequency of a SG wave. Note that the SG field 
resonances seen in Fig. ^ obey, according to condi- 
tion J4.10D , the relation uj/uj p < l/y/2 ~ 0.71. When 
the radius of the microsphere becomes substantially 



smaller than the wavelength Xp, then distinct peaks 
are only seen for a few lowest-order resonances [Rup- 



pin (1982D ; |Agarwal et al. (1983|) [. It is worth not- 



ing that, in contrast to dielectric matter, a large 
absorption in metals can substantially enhance the 
near-surface divergence of the decay rate, Eq. ( |5.8[ ), 
which is in agreement with experimental observa- 
tions of the fluorescence from dye molecules near a 
planar metal surface | prexhage (19 74| ) ] . 

8 QUANTUM CORRELATIONS 

Let us finally address the problem of spontaneous 
decay of an atom in the case when there is a second 
atom that is resonantly dipole-dipole coupled to the 
first one. Similarly to single-atom spontaneous de- 
cay, the dipole-dipole interaction can be controlled 
by the presence of macroscopic bodies. Various as- 
pects of the problem have been discussed for bulk 
material [ Knoester et al. (1989| )1, photonic crystals 
[[Kurizki et al. ( 1988j ); [Kurizki (1990|) ; [lohn et al. 



T995[ ); |Bay et al. (1997]); |Rupasov et al. (19970] 
optical lattices [ Goldstein et al. (1997b| ); |Guzman el 



1998|)] , planar ca vities |Kobayashi et al. (1995a ); 



garwal et al. (1998] )[ (and unspecified cavities [Kur 



|izki et al. (1996|); poldstein et al. (1997a )]), and mi- 
crospheres [ Agarwal et al. (2000 )]. In particular, res- 
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onant energy transfer realized through dipole-dipole 
interaction has been studied theoretically for bulk 
material [ Juzeliunas et al. (1994a]Jh|)], mi crospheres 
flGersten et al. (1984Q ; |Druger et al (19871) ; |Leung ej 



sphere (Appendix 



al. 
al. 



(1988[) ], and planar microcavities [Kobayashi e 



1995a, l])], and experimentally for droplets | Fola 



et c l. (1985 )1 an d planar microstructures [ Hopmeier] 



et 4/. (1999| ); |Andrew et al. (200Cp ], with potential 



for enhanced photon-harvesting systems and optical 
networks. 

Interatomic interaction can give rise to nonclassi- 
cal correlation and may be used for entangled-state 
preparation, which has been of increasing interest in 
the study of fundamental issues of quantum mechan- 
ics and with regard to application in quantum infor- 
mation processing. The entanglement, which is very 
weak in free space, may be expected to be enhanced 
significantly in resonator-like equipments. Propos- 
als have been made for entangling spatially sepa- 
rated atoms in Jaynes-Cummings systems through 
sequential or simultaneous strong atom-field cou- 



pling Kudryavtsev et al. (1993 


); Phoenix et al. 


(19§); 


Cirac et al. (1994 


); Freyberger et al. (1996|); 


Perry (1996|); picnio et al. (1999 


); Beige et al.\ 



Unfortunately, the Jaynes-Cummings model does 
not give any indication of the influence on entan- 
glement of the actual properties (such as form and 
intrinsic dispersion and absorption) of the bodies in 
a really used scheme. From Section [3.3| we know that 
not only the evolution of a single atom but also the 
mutual evolution of two atoms is fully determined 
by the Green tensor according to the chosen config- 
uration of macroscopic bodies. The formalism thus 
renders it possible to examine interatom quantum 
correlations established in the presence of arbitrary 
macroscopic body. 



8.1 Entangled-state preparation 



Weak Coupling Let us consider, e.g., two atoms 
near a microsphere of the type studied in Sec- 
tion ^. To be more specific, let as assume that 
the (two-level) atoms are of the same kind, that 
they are located at diametrically opposite positions 
(r J 4 = — r^), and that their transition dipole mo- 
ments are radially oriented. 

Obviously, the conditions ( |3.44 ) and ( 3.45| ) are ful- 
filled for such a system, so that from Eqs. (3.51) and 
( 3.53| ) together with the Green tensor for a micro- 



F: 



one then finds that 
™ p(l + !)(« + !) 

[ji(k A r A )+l3^(u JA )h\ 1) (k A r A ) 



xh?\k A r A ) [l=F(-l/ 



When atom A is initially in the upper state and atom 
B is accordingly in the lower state, then the two su- 
perposition states |+) and |— ), Eq. (|3.49j ), are ini- 
tially equally excited [C+(0) =CL(0) =2 - s]. If the 
atomic transition frequency coincides with a micro- 
sphere resonance, the most significant contribution 
to the single-atom decay rate T- 1 , Eq. (|7.1|) , comes 
(for sufficiently small atom-surface distance) from 
the corresponding term in the /-sum, i.e., 

r x ~ |r /(/ + i)(2/ + i) 



x Re 



h!t\k A r A ) 
k A r A 



1 2 



and Eq. ([S.ip can be approximated as follows: 



1T(-1) 



(8.3) 



Hence T_ S> T + (T + 3> r_) if I is even (odd), i.e., the 
state |— ) (|+)) decays much faster than the state |+) 

(I-))- 

Consequently, there exists a time window, dur- 
ing which the overall system is prepared in an en- 
tangled state that is a superposition of the state 
with the atoms being in the state |+) (|— )) and the 
medium-assisted field being in the ground state, and 
all the states with the atoms being in the lower state 
\L) and the medium-assisted field being in a single- 
quantum Fock state. The window is opened when 
the state |— ) (|+)) has already decayed while the 
state \L) emerges, and it is closed roughly after the 
lifetime of the state |+) (|— )). As a result, the two 
atoms are also entangled to each other. The state is 
a statistical mixture, the density operator of which 
is obtained from the density operator of the over- 
all system by taking the trace with respect to the 



medium-assisted field. Within approximation (8.3) 
it takes the form of 

p A ~ |C ± (t)| 2 |±)<±| + [1 - \C ± {t)\ 2 ] \L)(L\, (8.4) 

where 



C±(t)^2- 1 / 2 e- r ± t . 



(8.5) 
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Figure 14: The dependence of the decay rates 
r + (solid line) and T_ (dashed line) on (a) the 
transition frequency and (b) the distance of the 
atoms from a microsphere is shown for two atoms 
at (with respect to a sphere) diametrically oppo- 
site positions, radially oriented transition dipole 
moments, and a single-resonance Drude-Lorentz- 
type dielectric [R = 10 \t', cop = 0.5 cot; 7 = 10~ 6 o>t; 
Ar B = Ar A >W- 2 X T ; (a) Ar A = 0.02A T ; (b) 
to A — 1-0501 to?]- 



Applying the criterion suggested by Peres (1996| ) 



it is not difficult to prove that the state (8.4) is in- 
deed inseparable. It is worth noting that the atoms 
become entangled within the weak-coupling regime, 
starting from the state \U a) (or \Ub)) and the vac- 
uum field. In the language of (Markovian) damping 
theory one would probably say that the two atoms 
are coupled to the same dissipative system, which 
gives rise to the quantum coherence. 

The frequency dependence of T± as given by 
Eq. 



is illustrated in Fig. 14(a) for a frequency 
interval inside a band gap, and the dependence on 
the atom-surface distance is illustrated in Fig. |i~4|(b). 
We see that the values of T + and T_ can be substan- 
tially different from each other before they tend to 
the free-space rate Tq as the distance from the sphere 



becomes sufficiently large. In particular, the decay 
of one of the states |+) or |— ) can strongly be sup- 
pressed [see the minimum value of T+ in Fig. ||(b)] 
at the expense of the other one, which rapidly de- 
cays. Note that T + also differs from T_ for two 
atoms in free space flDeVoe et al. (19961 )]. However, 



the difference that occurs by mediation of the micro- 
sphere is much larger. For example, at the distance 
for which in Fig. 14(b) T + attains the minimum 
the ratio T_/r + ~ 67000 is observed, which is to be 
compared with the free-space ratio T_/r + ~ 1.0005. 
The effect may become even more pronounced for 
larger microsphere sizes and lower material absorp- 
tion, i.e., sharper microsphere resonances. Needless 
to say that it is not only observed for SG waves con- 



sidered in Fig. 14, but also for WG waves. 



Strong Coupling Entangled-state preparation in 
the weak-coupling regime has the advantage that it 
could routinely be achieved experimentally. How- 
ever, the value of |C+(t)| 2 in Eq. ( |3,4D (or the value of 
|C_(t)| 2 in the corresponding equation for the state 
|— )) is always less than 1/2. In order to achieve a 
higher degree of entanglement, the strong-coupling 
regime is required. 

Let us assume that the two atoms are initially in 
the ground state and the medium-assisted field is 
excited. The field excitation can be achieved, for ex- 
ample, by coupling an excited atom D to the micro- 
sphere and then making sure that the atomic excita- 
tion is transferred to the field (cf . Section |3.l| ) . If the 
atom D strongly interacts with the field, the excita- 
tion transfer can be controlled by adjusting the inter- 
action time. Another possibility would be measur- 
ing the state populations and discarding the events 
where the atom is found in the upper state. Here 
we restrict our attention to the first method and as- 
sume that all three atoms D, A, and B strongly 
interact with the same microsphere resonance (of 
mid- frequency coc an d line width Acoq)- According 
to Eq. ( 3.26| ), the upper-state probability amplitude 
Cjj d (t) of atom D reads 

C UD (t) = e -A"c(t+A*)/2 cos [ nD (t + At y 2 ], (8.6) 



with being given according to Eq. ( |3.27| ). For 

At = ir/n D , (8.7) 

the initially (i.e., at time t = — At) excited atom D 
is at time t = in the lower state [Cu D (0) = 0]. 

From the preceding subsection we know that 
when the resonance angular-momentum number I 
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is odd (even), then the state |+) (|— )) "feels" a 
sharply peaked high density of medium-assisted field 
states, so that a strong-coupling approximation of 
the type ( 3.24 ) applies. The state |— ) (|+)), in 
contrast, "feels" a flat one and the (weak-coupling) 
Markov approximation applies. Assuming atom A 
(or B) is at the same position as was atom D, from 
Eqs. ( pT47D , (pD, and Q) we then find that 



c±(t) 



-Auj c (t+ir/n D )/2 



sin(fi±t/2) 



[with Q± according to Eq. ( |3.27| )] and 

C T (t)~0 (8.9) 

(the sign of C_ (t) in Eq. ( |§lj| ) is reversed if atom B is 
at the same position as was atom D). Note that Vt± 
= 2 2 ft, because of Eq. fl8.3|) . The two-atom entan- 
gled state is again of the form given in Eq. ( p.4[ ) , but 
now the weight of the state |+) (|— )) can reach val- 
ues larger 1/2, provided that the resonance linewidth 
Auic is small enough. 




8.2 Violation of Bell's inequality 

Violations of Bell's inequalities provide support to 
quantum mechanics versus local (hidden-variables) 
theories [[Bell (1965| ); |Clauser et al. (1969| )1. De- 
spite outstanding progress^ in the test of Bell's in- 
equalities, a decisive experiment to rule out any lo- 
cal realistic theory is yet to be performed JVaidman 
(2001 )1, and the problem continues to attract much 



attention. Though entangled states of spatially sep- 



arated atoms in a cavity have been observed [ Hagley 
et al. (1997)], a test of Bell's inequalities for such a 



system has yet to be realized. 

The Bell's inequality for spin systems can be writ- 
ten in the form of [Bell (1965| ); |Clauser et al. (1969| )1 



B S = \E(9 1 ,9 2 )-E(9 1 ,9' 2 ) 
+ E(9' 1 ,9 2 )+E(9[,9' 2 ) 



< 2, 



where 



E(e u e 2 



cos9a A + sin a y A . 



(8.10) 

(8.11) 
(8.12) 



When the atomic state \ua,ub) is not populated, as 
it is the case for a state of the type ( |8.4| ), it is not 
difficult to prove that 



E(e 1 ,e 2 ) = E(9 1 -e 2 ,o). 



5.13) 



Figure 15: The dependence on time of B$ is 
shown for two atoms at (with respect to a mi- 
crosphere) diametrically opposite positions, radially 
oriented transition dipole moments, and a single- 
resonance Drude-Lorentz-type dielectric [R = 10 At; 
ujp = 0.5 ujt', Are = Ar A = 0.02 At; oj a = 1-0501 lot', 
r = 10- 6 u>r; M D = ft; j/oj t = W- 6 (solid line), 
10 -5 (dashed line)], (b) ft/Awe versus Ar A for 
7/wt = 10 -6 (Ar A > 10~ 3 A T ). The inset shows the 
variation of the first maximum value of B$ in (a) . 



Let us choose 



C — U\ — U 2 — U 2 — — (7j — U 2 . 

The inequality ( |8.10| ) thus simplifies to 
B s = \3E(0,0) -£(30,0)1 < 2. 



.14) 



1.15) 



An entangled state of the type (8.4) can only 
give rise to a violation of the Bell's inequality if 
|C+(*)| 2 > —0.71 jBeige et al. (200CQ 1, which 
cannot be achieved in the weak-coupling regime, 
Eq. (8.5). It can be achieved, in contrast, in the 



5 F or recent experiments us ing; photons, see Weihs et al 
f!99^); iKuzmich et al. (2000| ), and using trapped ions, see 
Rowe et al. (200l[). 
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strong-coupling regime, Eq. 
E(6,0) = cos 9 \C±(t)\ 2 



where 



cos e -A«c(*+Vn D ) sin 2 (nt/V2) . (8.16) 



Substitution of this expression into in Eq. (|0|) 
yields, on choosing # = 7r/4, 

B s = 2V2e- AuJc( - t+n / nD hin 2 (nt/V2), (8.17) 

which clearly shows that Bs > 2 becomes possible as 
long as Auj c(t + tt/Qd) *C 1. Examples of the tem- 
poral evolution of Bs for two atoms near a dielectric 
microsphere are shown in Fig. |D^(a). In Figure 15(b) 
the dependence of the ratio 0,/Auc ° n the distance 
of the atoms from the sphere is plotted. The strong- 
coupling regime can be observed for distances for 
which 0/Ao;c , S>l is valid. The inset reveals that 
the maximum value of Bs decreases with increasing 
atom-surface distance and reduces below the thresh- 
old value of 2 still in the strong-coupling regime. 

9 SUMMARY 

We have studied spontaneous decay in the pres- 
ence of dispersing and absorbing macroscopic bod- 
ies, basing on quantization of the (macroscopic) elec- 
tromagnetic field in arbitrary linear, causal media. 
The formalism covers both weak and strong cou- 
plings and enables one to include the material ab- 
sorption and dispersion in a consistent way, without 
restriction to a particular frequency domain. It re- 
places the standard concept of orthogonal-mode de- 
composition, which requires real permittivities and 
thus does not allow for material absorption, with a 
source-quantity representation in terms of the classi- 
cal Green tensor and appropriately chosen bosonic- 
field variables. All relevant information about the 
bodies such as form and intrinsic dispersion and ab- 
sorption properties are contained in the Green ten- 
sor. 

The formalism has been applied to study spon- 
taneous decay of a single atom in the presence of 
various absorbing and dispersing macroscopic bod- 
ies, including open configurations such as bulk and 
planar half space media, and closed configurations 
such as a spherical cavity or a microsphere. Ab- 
sorption can noticeably influence spontaneous decay. 
So, the decay rate in absorbing bulk material takes 
a much more complicated form than one would ex- 
pect from the simple product form that is commonly 
used for nonabsorbing matter. The decay rate of an 



atom located very near a planar surface shows that 
due to material absorption the decay rate drastically 
rises as the atom approaches the surface of the body, 
because of near-field assisted (nonradiative) energy 
transfer from the atom to the medium. In fact, this 
is valid for an atom that is sufficiently near an arbi- 
trary body, because for short enough atom-surface 
distances, any curved surface can be approximated 
by a planar one. 

Spontaneous decay can strongly be influenced by 
field resonances that can appear due to the presence 
of macroscopic bodies, depending on whether the 
atomic transition frequency is tuned to a field reso- 
nance or not. In particular, the decay process can 
be mainly radiative or nonradiative, which depends 
on whether the radiative losses due to input-output 
coupling or the losses due to material absorption 
dominate. In particular, to understand what hap- 
pens when the atomic transition frequency is inside 
a band-gap zone of a body, inclusion in the study of 
material absorption is necessary. 

Finally, spontaneous decay of (two) dipole-dipole 
coupled atoms in the presence of macroscopic bodies 
offers the possibility of entangled state preparation 
and verification of the violation of Bell's inequalities. 
Whereas entangled states can already be prepared 
in the weak-coupling regime, violation of Bell's in- 
equalities requires the strong-coupling regime, the 
ultimate limits being given by material absorption. 
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A THE GREEN TENSOR 

A.l Bulk Medium 

For bulk material the Green tensor reads as 



[p=r-r 



G(ry,u J )=[V^^ + I q ^)]^ m -,(A.l) 



where 



= \J s{oj) uj/c. 



(A.2) 
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It can be decomposed into a longitudinal and a 
transverse part, 

G(r,r',u) =Gll(r,iV) + G ± (r,r / ,cu), (A.3) 

where 



G"(r,r',u/) 



1 



47TQ' 2 



47T 



3p<g>p\ 1 



and 



G ± (r,r',o;) 



1 



1 



4-7rg 2 



3p<g>p\ 1 



(A.4) 



1 



I 



1 



3i 3 \ p8p 



K qp (qp) 2 (qp) 3 J p 2 
In particular, from Eq. ( |A.5| ) it follows that 
ImG 1 (r,r,u) = lim Im G (r, r', lj) 

r'— >r 



(A.5) 



67TC 



n R (u)I. 



(A.6) 



A. 2 Spherical Multilayers 



The Green tensor of a spherical structure, consisting 
of N concentric layers, enumerated from outward in 
(the outmost layer labeled as layer 1, the innermost 
layer as layer AT), can be decomposed into two parts 

G(r, r', u) = G (s) (r, r', Lo)5 fs + cW(r, r', u), 

(A.7) 

where G^ s '(r,r' ,u) represents the contribution of 
the direct waves from the source in an unbounded 
space, and G^ s \r,r' ,uj) is the scattering part that 
describes the contribution of the multiple reflection 
(/ = s) and transmission (/ / s) due to the presence 
of the surfaces of discontinuity (/ and s, respectively, 
refer to the regions where are the field and source 
points r and r'). In Eq. QA.7 ), G^ s ' is nothing but 
the bulk- material Green tensor ( A.l ). In the (local) 
spherical coordinate systems it reads as [see, e.g., [Li| 



et al. (1994|) 1 



G W( 



r,r ,u) 



S(r 



,i_hs^ST^ [ (l-m)\ 
+ 4vr 2^2^ 2^1 m+i) Q+m)\ [ ' 

e 1=1 m=0 k v ' v y 



(i) 



+ N^(r,y®N eim (r',fc s 



(A. 



if r>r', and G^ s \r,r',u) = G^(r',r,u>) if r<r'. 
The scattering partG^(r, r',u) in Eq. (|A.7[) reads 



[ |Li ef g (19941) 1 
G^(r,r', 



oo i 



47T 



Z=l m=0 
r(l) 



1(1 + 1) (l+m)\ 



(1 ~ ^Af)Mej ra (r, kf) <8> [M«j m (r', fc s ) 
x (1 - + M^Jr'^Xl - 5^)13? 

o 

+ (1 - 5 /A r)N ( ^ m (r, fc/) ® [Ne |m (r', fc a ) 

o ° 

x (1 - 5 sl )A?(u) + N^Jr', fc s )(l - <W)S/»] 

o 

+ (1 - 6 fl )Me lm (r, k f ) [Me^(r', k s ) 

J o J L o 

x (1 - 5 sl )C/V) + M^Jr', fc.)(l - <W)2?f M] 

o 

+ (l-<y/ 1 )Ne, m (r,A;/)® [N. (m (r', fc s ) 
x (1 - S^cFlu) 

+ Ni 1 l i (r',k s )(l-5 sA f)VF(u;)])), 



(A.9) 



where 



(A.10) 



and M and N represent TM- and TE-waves, respec- 
tively, 



Me nm (r, k) = T ^-j n (kr)P™(cos 9) ( Sm ) (W>)e e 
° sm # \cos / 

. . dP™(cos9) /cos\ , ,. . . ^, 

a9 \ sm / 



t (r, fc) = n \ + 1 ^ n (fcr)P^(cos (9) ( C ° S I ! /m > )e,. 
fcr V sm , 



1 d[rj n (£r)] 



^—E P n(cOs9)[ 

sm 9 \cos 



dP™( cos 9) /cos N 
aft 1 \s 



, sm 



(m(f))e 4 



(A.12) 



with j n (x) being the spherical Bessel function of the 
first kind and P™ l (x) being the associated Legendre 
function. The superscript (1) in Eq. ( |A.9| ) indicates 
that in Eqs. ( |A.11| ) and ( |A.12|) the spherical Bessel 
function j n (x) has to be replaced by the first- type 
spherical Hankel function hn\x). 
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are to be found from the coupled recurrence equa- 
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[with / and s being taken according to Eq. 
Here, the coemcients are redefined as 
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with 



Jil = 3n(kiRl), 

H a = h^(kiRi), 

T , 1 d[P3n(p)] 
J u 



P dp 

ld[ph { n\p)} 



p=kiRi 



P dp 



(A.24) 
(A.25) 
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p=kiRi 



The coefficients in the first and the last layers can 
be found immediately from Eqs. ( |A.13 ) and ( [A. 14 ). 
The rest can be obtained by again using recurrence 



equations (A.13) 



A. 2.1 Two-layered medium 

For a sphere of radius R (including the special cases 
of an empty sphere in an otherwise homogeneous 
medium and a material sphere in vacuum) we have 
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(A.29) 

(A.30) 
(A.31) 



A. 2. 2 Three-layered medium 

For three-layered media of radii i?i and i?2 (i?i > 
R2), the spherical cavity presented in Fig. || in par- 
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ticular, we have that |Li et al. (1994|)] 



e n=l i=0 



n(n+l) (n+Z)! 
^MMlJ, (r J fci)®M« ni (r / ,fc 3 
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